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We have evaluated the use of phase synchronization to identify resting state functional 
connectivity (RSFC) in the language system in infants using functional near infrared 
spectroscopy (fNIRS). We used joint probability distribution of phase between fNIRS 
channels with a seed channel in the language area to estimate phase relations and to 
identify the language system network. Our results indicate the feasibility of this method 
in identifying the language system. The connectivity maps are consistent with anatomical 
cortical connections and are also comparable to those obtained from functional magnetic 
resonance imaging (fMRI) functional connectivity studies. The results also indicate left 
hemisphere lateralization of the language network. 

Keywords: functional near-infrared spectroscopy, resting state functional connectivity, phase synchronization, 
language network, BOLD signal 



1. INTRODUCTION 

Functional Near Infrared Spectroscopy (fNIRS) is an optical 
method for brain functional imaging. Using near infrared light 
with high penetration depth in human tissue, fNIRS can measure 
the neuronal activity in response to a task or an environmen- 
tal stimulation through neurovascular coupling in superficial 
areas of the brain. fNIRS has been used for functional stud- 
ies as a portable, less expensive and less restraining alternative 
to functional Magnetic Resonance Imaging (fMRI) in different 
task based brain functional studies. One of the more recent areas 
of interest in both fNIRS and fMRI is the study of interaction 
between different cortical areas through their intrinsic neuronal 
signaling (Biswal et al, 1995; White et al, 2009). This intrin- 
sic signaling appears in the form of slow varying spontaneous 
fluctuations in the blood oxygen-level dependent (BOLD) signal 
in the absence of stimulation. These fluctuations are correlated 
between brain areas that are anatomically and functionally con- 
nected and have been used to map brain functional networks 
such as sensorimotor, visual, and auditory as well as higher net- 
works such as language and attention (Beckmann et al., 2005; 
Zhang and Raichle, 2010). These networks are often mapped from 
data collected at rest and are therefore referred to as Resting State 
Functional Connectivity (RSFC) maps. 

The RSFC analysis can also potentially identify changes in 
intrinsic neural activity as a result of disease in some neurological 
and psychiatric conditions. Changes in connectivity strength in 
different brain networks have been observed in conditions such 
as autism (Cherkassky et al, 2006), depression (Anand et al., 



2005), Alzheimer disease (Li et al, 2002), and attention-deficit 
hyperactivity disorder (Tian et al., 2006). 

Given the advantages of fNIRS, different brain networks have 
been investigated through RSFC using fNIRS. One of the most 
common methods for analyzing brain network connectivities 
using RSFC in fNIRS is cross correlation (Zhang and Raichle, 

2010) . In cross correlation, an fNIRS channel is selected as the 
seed channel and the correlations of the signal in all other chan- 
nels with the seed channel are calculated. The objective is to 
find the cortical areas whose resting state fluctuations are simi- 
lar to that of the seed channel. Cross correlation based functional 
connectivity has been used in conjunction with fNIRS to derive 
connectivity maps in different brain networks (White et al., 2009; 
Lu et al., 2010; Duan et al., 2012). Cross correlation and cross 
coherence of changes in fNIRS channels in particular have been 
used on infants to study functional connectivity. In a study by 
Homae et al, the functional connectivity in 3 month old infants 
with and without audio stimuli was investigated (Homae et al., 

2011) . The study showed that high correlation exists between 
neighboring channels and also their corresponding channels in 
the contralateral hemisphere during resting state. One drawback 
of correlation based connectivity is that it can be sensitive to 
detection of spurious connection as a result of presence of cross 
talk between channels, systemic interference or noise (Tass et al., 
1998). 

The phase of the oxygenation changes in fNIRS has been 
shown to contain information on functional connectivity. In a 
study by Taga et al., strong phase synchronization was observed 
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between cortical regions during sensory processing in young 
infants (Taga et al., 2011). The phase information was also used 
in this study to determine the timing and flow of activity in 
the developing brain. The phase relationship between oxygenated 
and deoxygenated hemoglobin in resting state in infants was 
investigated in another study by Taga et al. (2000). Spatially syn- 
chronized oscillations were also observed in this study in the mea- 
surement channels over the occipital cortex (Taga et al, 2000). 
Phase based methods have also been used in other modalities such 
as EEG and MEG for functional analysis (Stam et al., 2007). 

Previous studies have investigated the neural connectivity in 
infants using fNIRS. However, there have been limited studies on 
the neural systems involved in speech processing in the newborn 
with some level of inconsistencies in the data regarding the later- 
alization of the network. In this paper, we have investigated the 
application of a new method based on the phase relation between 
spontaneous oxygenated hemoglobin changes in fNIRS channels 
and have used it as a measure of functional connectivity to iden- 
tify the speech processing neural systems in newborn infants. 
Compared to methods based on signal amplitude such as cross 
correlation, phase is much less sensitive to noise and interferences. 
It also does not require the assumption of stationarity for the sig- 
nals. The phase synchronization is not equivalent to coherence or 
frequency synchronization and is an independent characteristic of 
the interrelationship between two processes (Tass et al., 1998). To 
evaluate the feasibility of this analysis method for detecting func- 
tional connectivity, we applied it to a study of processing of speech 
vs. non-speech in newborn human infants. This type of compari- 
son is of particular value for the question being asked as there is a 
growing body of evidence on the brain areas involved in language 
processing in neonates, but less on the underlying connectivity. 

2. MATERIALS AND METHODS 
2.1. fNIRS DATA 

The fNIRS data were collected from newborn infants at BC 
Children's Hospital, Vancouver, Canada, during a separate lan- 
guage perception study (May et al, 2012). Informed consent 
was obtained from parents when the experiment was being con- 
ducted. The study design was approved by the ethics committee 
of the University of British Columbia. The experiment design and 
setup are shown in Figure 1. A total of 19 subjects were used in the 
analysis (mean age 1.4 days, 7 females, age range 0-4 days) out of 
which two subjects were excluded due to severe artifacts in the 
signals and poor data quality resulting from optode displacement 
during data collection. During the experiment, audio stimuli were 
administered to subjects while the subjects were in state of quiet 
rest or sleep. The audio stimuli consisted of blocks of sentences 
in Spanish and Silbo-Gomero. Silbo-Gomero is a whistled lan- 
guage that is a surrogate language of Spanish (Carreiras et al., 
2005; Rialland et al., 2005) It uses whistles rather than speech, 
and was developed by shepherds in the Canary Islands to com- 
municate across long distances. Spanish and Silbo-Gomero were 
selected as both are unfamiliar to the infants, while one is a spoken 
language and the other is not. Each block was 15 s long followed 
by 25-35 s of silence. A total of 8 blocks for each stimulus were 
presented in which each block consisted of continuous speech. 
The total experiment time was 22-25 min. 



8 Blocks 



I Silb 1 1s 



Silence Silb 2 I Silence Silb 3 Silence 



J Silb 4 [ 



Silence 



15 s 25 s 15 s 35 s 15 s 35 s 15 s 25 s 




ETG-4000 
NIRS 
System 



Markers 



t 



FIGURE 1 | fNIRS experiment setup. 



The subjects' brain hemodynamic response was monitored by 
a 24 channel fNIRS device (Hitachi ETG-4000 machine with 695 
and 830 nm lasers at a power of 0.75 mW, interoptode distance of 
3 cm and sampling rate of 10 Hz). Two chevron shaped optode 
holders were used for securing nine 1 mm fibers to the head. 
There were a total of 4 detector and 5 source fibers on each holder 
resulting in 12 recording channels per holder. Figure 2 shows the 
placement of optodes on the subject's head. Surface landmarks 
(ears or vertex) were used for the placement of the probe holder 
over the infants' perisylvian area of the scalp. Channels 1 1 and 12 
in LH and 23 and 24 were ideally placed above the infant's ear. A 
stretchy cap was used to secure the holders on the infants' head. 

2.2. DATA ANALYSIS 

In order to determine the phase relation between channels, we 
first extracted the phase of the signal in each channel using the 
Hilbert transform. The Hilbert transform converts a real valued 
signal to a complex one, known as analytic signal, whose real part 
and phase correspond to the original signal and its derived phase, 
respectively (Oppenheim et al., 1989). The Hilbert transform of 
signal x[n] in the frequency domain is defined as (Oppenheim 
etal., 1989): 

Y(<J M ) = -jsgn( ti >)(x(<J< a )) (1) 

where X (e 7 ™) is the Fourier transform of x[n] and sgn(a>) is the 
sign function having value of 1 for 00 > 0 and -1 for co > 0. The 
analytic signal can then be written as 



x a [n] = x[n] + jy[n] 



(2) 



where y[n] is the inverse Fourier transform of Y(e JW ). 

We used the joint probability distribution of the phases across 
channels to describe their connectivity. A common model for 
probability distribution of phase which is the circular analog of 
the Gaussian distribution is the Von Mises distribution. The Von 
Mises probability density function (pdf ) is defined as (Mardia and 
Jupp, 2000): 



/(0IM-.K) 



1 



k cos(8 — ^l) 
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(3) 
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FIGURE 2 | Optode placement on the head. Blue squares and red dots indicate detectors and sources, respectively while the dashed lines indicate the 
measurement channels and the numbers indicate the channel number. 



where 6 is an angle defined in the interval [— jt, Tt) and Iq(k) is 
the modified Bessel function of order 0. The parameter k is the 
equivalent of the covariance for the Gaussian distribution and u, is 
the expected value of the angle. The probability density function 
of the signal phase in channel m conditioned on that of channel n 
can therefore be written as: 



/(9 m - 0„|y,, k) 



1 



jc mlt cos(9 m — 9 n — 



(4) 



K m „ describes the intensity of the phase correlation between 
signals in channels m and n. In other words, it shows how much 
the prior information of 6„ affects distribution of Q m . The first 
moment of the distribution given in Equation 3 can be calculated 

as (Mardia and Jupp, 2000) 



mi =E[ei e ] 



f e> e f(Q\[L,\c)dQ 
J— it 



(5) 



Using the first moment, one can estimate parameter k by numer- 
ically solving the optimization problem 
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(7) 



where N is the number of samples in the data segment and (j), is 
the measured phase of the signal at time point i. Parameter |x can 
then be estimated using 



N ^ 



(8) 



A close relationship exists between parameters of the distri- 
bution and the phase locking value (PLV) which is a common 
measure used in EEG signal processing to detect functional 
connectivity through synchronization between channels. PLV is 
related to the distribution through the magnitude of the first cir- 
cular moment of the phase distribution (Canolty et al, 2012) 



PLV„ 



;|> m -e„)]| 



ft 00 
JbOO 



(9) 



One advantage of using Von Mises distribution over phase 
locking for connectivity analysis is that once the parameters are 
estimated, one would have the distribution function and can re- 
sample from the distribution to determine the significance levels. 
Also, the preferred phase difference is not available in PLV. 

To evaluate fNIRS functional connectivity, we first calculated 
the phase from all fNIRS channels using the Hilbert transform 
as described earlier. Ideally, we would be interested in phase rela- 
tions between channels when subjects are not exposed to any type 
of stimulation to reveal intrinsic network activities. It was shown 
in an earlier study, however, that the infant brain produces no 
significant response in the language network to Silbo-Gomero 
stimuli (May et al., 2012). We therefore used the fNIRS data dur- 
ing Silbo-Gomero stimulation as an alternative to resting state. 
The BOLD spontaneous fluctuations are concentrated at frequen- 
cies of less than 0.1 Hz (Biswal et al., 1995). Therefore, the signals 
were first filtered with an infinite impulse response bandpass filter 
(IIR) between 0.02 and 0.08 Hz to extract spontaneous hemo- 
dynamic activities and reject other interferences. This frequency 
range is comparable to those used in other studies investigating 
RSFC using fMRI and fNIRS (Biswal et al, 1995; Duan et al, 
2012). 

fNIRS data in general can be contaminated with motion arti- 
facts as the result of subjects' spontaneous movements. These 
artifacts create interference in the form of highly correlated phase 
changes in fNIRS channels, especially in spatially close chan- 
nels. This interference results in a very high phase correlation 
and can obscure underlying phase connections between channels. 
Even though filtering of the motion artifacts is possible, in order 
to minimize possibility of introducing any inter-dependence 
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between channels, no artifact removal procedure was applied. 
Instead, the channels for all subjects were inspected visually and 
artifact contaminated regions within the Silbo-Gomero stimula- 
tion window were marked. An artifact free segment of the data in 
each channel was then selected for the analysis and the phase of 
the selected signal segments were then derived using the Hilbert 
transform. Since the brain shows no response to this stimula- 
tion type, the stimulation onsets were ignored and the segments 
were selected independent of the stimulation onsets. The seg- 
ments contained variable number of stimulation blocks and their 
length ranged from 50 to 220 s. 

The channel with the highest activation in the grand average 
for the Spanish stimulation task during the original study in the 
left hemisphere was selected as the seed channel for the RSFC 
analysis. The joint phase distribution of the seed channel and all 
other channels was then estimated by calculating the phase dif- 
ference between the seed channel and other channels and then 
estimating k,„„ and |x m „ using Equations 7 and 8. We used a sim- 
plex derivative-free method to solve Equation 7 and derive \c mn 
(Lagarias et al., 1998). The analysis was performed in MATLAB 
(Mathworks MA, USA) and the phase coupling estimation tool- 
box developed by Cadieu et al. was used for parts of the analysis 
(Cadieu and Koepsell, 2010) 1 . 

The analysis was performed on oxygenated hemoglobin 
(HbCh) changes only. Previous studies on the application of 
fNIRS to detect language network activity and connectivity have 
shown that the oxygenated hemoglobin is more sensitive to 
regional cerebral blood flow changes than the deoxygenated 
hemoglobin with the equipment used here (Hoshi, 2007; Gervain 
et al, 2008; Lu et al, 2010). 

To examine the validity and reliability of the connectivity 
information derived with this method, we divided the subjects 
randomly into two groups and evaluated the connections for each 
group, similar to the approach used in Zhang et al. (2010). We 
then compared the correlation of the connectivities between the 
groups. 

Other studies have suggested that the language network is left 
lateralized. We verified this in the network derived using our 
method. The lateralization was quantified using (Binder et al., 
1996; Zhang et al., 2010) 



LI ■■ 



1 



M/2 



E 



= 1 K L 



+ 4 



(10) 



where M is the total number of channels, k' l is the value of k; 5 
in which i is the channel number and s is the seed channel in 
the left hemisphere, is the value of the same parameter with 
the channel symmetric to i in the right hemisphere. The signifi- 
cance level of the calculated lateralization index is then evaluated. 
The lateralization index results in a number between —1 and 
1 with more positive numbers indicating higher degrees of left 
lateralization. 

As the final step, we defined 4 regions of interest (ROI), 2 inside 
the language network and 2 outside the network and evaluated the 

1 http://redwood.berkeley.edu/klab/pce.html. 



connection strengths in these areas. In particular, channels 6 and 
7 were selected inside the network, based on our prior knowledge 
that the physical area they cover is in the language network, and 
channels 1 and 12 outside the network with the optode configura- 
tion used in the current study. Channel 1 is over the frontal areas 
while channel 12 covers the temporal area. The choice of these 
channels as being outside the language area is justified by the fact 
that they showed no significant activation in response to native 
language or Spanish (May et al., 2012). 

3. RESULTS 

Figure 3 shows the Hb02 signal from channels 7 and 9, the seed, 
for a typical subject. Qualitatively, the histogram of the phase 
at different time points for both channels does not show a clear 
dominant phase range as shown in Figures 3B,C. However, the 
joint distribution histogram shown in Figure 3D has a sharp peak 
focused around the mean phase difference. This is also evident in 
the Von Mises distribution function plot derived from estimated 
parameters in each case (shown in red). The estimated distribu- 
tion parameters are also indicated in the fig ure. In the case of the 
phase histogram for individual channels (Figures 3B,C), the val- 
ues of estimated K m „ are much smaller than in the conditional 
distribution (Figure 3D). This indicates a high phase relationship 
between the two channels and is interpreted as connectivity. 

Using the method described earlier, the group level resting 
state functional connectivity map with channel 9 chosen as the 
seed channel was derived and is shown in Figure 4. The detected 
network includes the areas known to be associated with lan- 
guage network including the superior temporal gyrus and Broca's 
area. The maps are also in agreement with those obtained for 
the language network in adults using correlation based fNIRS 
connectivity studies (Zhang et al., 2010). 

The connectivity maps resulting from the 2 random subgroups 
are shown in Figure 5. The maps for the two subgroups cover 
similar areas in both hemispheres. The correlation between indi- 
vidual connections in the two subgroups is shown in Figure 6 
(Pearson correlation r = 0.6). 

Figure 7 shows the results of ROI connectivity analysis where 
the connection strength between the seed channel and 2 chan- 
nels in the language area (6 and 7) is compared with that with 
two channels outside the language system (channels 1 and 12). 
Analysis of variance indicates significant difference between the 
connections (AN OVA p < 0.01) inside and outside the language 
network. In particular, channels 6 and 7 connections are not dif- 
ferent while they are both higher than that of channels 12 and 1 
in the temporal and frontal areas, respectively. Results of multiple 
comparison test are shown in Table 1 (Tukey's test). 

The results of the lateralization analysis are shown in Figure 8. 
The average lateralization index is 0.172 and is significantly 
different from zero (1 sample f-test, p < 0.001). Here, the lat- 
eralization index is also compared for all subjects between the 
language network and a control case. The control network is 
created by choosing channel 11 as the seed channel. There 
is a significant difference in lateralization index between the 
language and control network (paired f-test p < 0.01). These 
results suggest strong left lateralization in the detected language 
network. 



Frontiers in Human Neuroscience 



www.frontiersin.org 



January 2014 | Volume 7 | Article 921 | 4 



Molavi et al. 



Analyzing functional connectivity using NIRS 



A 0.4 




= -0.42 



Phase (rad) 1 2 3 3 phas ° ^ 1 1 phas 3 ^ 



FIGURE 3 1 Filtered Hb0 2 signal recorded from two channels channel 9] along with the joint phase distribution (D). The 

with high degree of connection (A) and the distribution of red curve in histograms shows the estimated probability 
the phases for each channel [(B) for channel 7 and (C) for distribution. 
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FIGURE 4 | Group level RSFC maps with channel 9 in the language area used as the seed channel. Left and right panels correspond to left and right 

hemispheres. 



4. DISCUSSION AND CONCLUSION between fNIRS channels with a seed channel in the language 

We evaluated use of phase synchronization to identify resting area to estimate phase relations and identify the language sys- 

state functional connectivity in the language system in infants tern network. Our results indicate the feasibility of this method 

using fNIRS. We used joint probability distribution of phase in identifying the language system. The connectivity maps are 
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consistent with anatomical cortical connections and are also com- 
parable to those obtained from fMRI functional connectivity 
studies (Perani et al., 2011; Mahmoudzadeh et al., 2013). The 
results indicate left hemisphere lateralization of the language 
network. 

The brain networks connectivity reveals information about 
underlying anatomical areas involved in a particular task. In 
some disease conditions, changes in cortical connections occur 
(Zhang and Raichle, 2010). Application of connectivity estimat- 
ing methods to fNIRS enables investigation of such changes 
in cases where use of fMRI is not possible, such as in infants 
and extends utilization of fNIRS in wider range of clinical 
applications. 



Use of fNIRS for analysis of functional connectivity offers sev- 
eral advantages over more traditional fMRI based connectivity 
analysis. Collection of fMRI data from infants and young children 
under resting condition can be challenging. In contrary, fNIRS 
is easily applicable to even newborns. Also, in cases where sub- 
jects to be tested are immobilized and can not be transferred to an 
MRI scanner, portable fNIRS systems can replace fMRI for con- 
nectivity analysis. One limitation compared to fMRI is the limited 
penetration depth which means connectivity analysis will be lim- 
ited to cerebral cortex. Use of phase information and in particular 
using the method described here has some advantages that can 
justify its use for functional analysis. In comparison with ampli- 
tude based methods such as cross correlation and cross coherence, 
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FIGURE 6 | Correlation between connection strengths in the two 
subgroups. 
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FIGURE 7 | ROI connections comparison between language area 
(represented by channels 7 and 6) and outside language area 
(represented by channels 1 and 12). The bars indicate the standard error 
of the mean. 



Table 1 | Pairwise comparison between selected ROIs inside and 


outside language network (Tukey's test). 




Channel Pair 


7-6 7-12 


7-1 12-1 


Mean Ak 


-0.15 1.65 


1.57 -0.09 


CI (95%) 


[-1.43 1.13] [0.37 2.94] 


[0.29 2.85] [-1.37 1.20] 



phase-based methods do not suffer from spurious connections 
as a result of common source. The method described here can 
be generalized to estimate the phase distribution function in a 
multivariate approach (Canolty et al., 2012). In this way, the con- 
nectivity parameter between all channels is estimated at once. 
This eliminates the limitation of bivariate methods in detection 
of spurious connections. 

One limitation of fNIRS is the lack of information of anatom- 
ical structure of participants' brains. Due to the lack of an 
established standard system for processing infant brains, some 
procedure to localize the channels to the brain may be neces- 
sary. The probe set configuration we used in this study primarily 
allows measurement of activation from anterior temporal to pos- 
terior temporal areas, but given infants' small heads, the anterior 
temporal areas likely extend into frontal regions and the posterior 
temporal likely extend into parietal regions. Nonetheless, given 
the variability among infant head sizes and the lack of agreed 
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FIGURE 8 | Lateralization index for the language network and 


comparison with control network. 







upon stereotactic maps of infant heads, we chose the most con- 
servative approach of using the general terms we have used, such 
as anterior and posterior temporal, to refer to the areas measured. 

Our results are comparable to similar studies in the literature. 
In particular Zhang et al. analyzed RSFC in the language system 
in adults using fNIRS (Zhang et al, 2010). Their results indicated 
significant RSFC between left inferior frontal and superior tem- 
poral cortices which are associated with language system (Zhang 
et al., 2010). Using fMRI, Fransson et al. studied resting state net- 
works in the infants brain (Fransson et al., 2007). They observed 
similar networks in the bilateral temporal/inferior parietal cor- 
tex which encompasses primary auditory cortex (Fransson et al., 
2007). In a study using fMRI on newborn infants, Perani et al. 
showed strong interhemispheric functional and structural con- 
nections between left and right temporal regions and also between 
the left and right frontal regions (Perani et al., 2011). Homae et 
al used the cross correlation and spectral coherence to investi- 
gate the functional connectivity in infants' brain (Homae et al., 
201 1). This study reports high correlation between adjacent chan- 
nels and the corresponding contralateral channels. This is in 
agreement with our findings of connectivity between neighboring 
channels in the language area in the left hemisphere and corre- 
sponding channels in the right hemisphere using the phase-based 
method described here. 

Other phase-based methods for analysis of functional connec- 
tivity have been reported earlier. In a study by Taga et al, the phase 
information of the oxygenated and deoxygenated hemoglobin 
was investigated in the presence of audio stimuli (Taga et al., 
2011). The relationship between the stimulation and the evoked 
responses were used to evaluate the temporal features of the 
response. One important difference of our work is the absence 
of the stimulation evoked changes for detection of functional 
connectivity. 

The presence of motion artifacts has a significant effect on the 
connectivity strength results using phase synchrony method pre- 
sented here. The motion artifact results in in-phase changes across 
affected channels which may result in stronger phase correlation 
compared to those resulting from spontaneous neuronal activity. 
Therefore, care must be taken to ensure segments being processed 
do not include motion artifacts. Various methods for removal 
of motion artifacts in fNIRS have been proposed (Cooper et al., 
2012; Brigadoi et al., 2014). In this particular study, we used visual 
inspection to remove motion contaminated data segments instead 
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of using an automated removal technique. An artifact removal 
method may introduce some level of distortion in the artifact- free 
regions of the signal and also may not fully remove the artifacts. 
Therefore, to avoid detection of connectivity as a result of such 
distortions, an automated approach was not employed. Saturated 
channels or channels that have lost coupling to tissue due to 
displacement will also have similar effect. 

Our fNIRS data was not collected during strict "resting" state. 
We used continuous blocks of data in which subjects were listen- 
ing to a non-speech audio stimulation. No significant activation 
compared to baseline was observed on this data and was there- 
fore used as the baseline. Some fMRI studies have used a similar 
approach for mapping RSFC. The study by Greicius et al. on 
default mode network in Alzheimer's disease patients for exam- 
ple, was performed during a low demand cognitive task (Greicius 
et al, 2004; Zhang and Raichle, 2010). An alternative approach 
in the case of task-related data is to regress out the task evoked 
response from the data before performing RSFC analysis (Fair 
etal.,2007; He etal., 2007). 

The fNIRS signal is known to contain systemic interferences. 
This includes interference from cardiac pulsation, respiration, 
cardiovascular autoregulation and heart rate variability. The 
frequency band for connectivity analysis must be chosen such 
that it includes the relevant variations caused by the neurovascu- 
lar coupling while rejecting the frequency bands containing these 
interferences. The cardiac interference in our study is around 2 Hz 
and the very low frequency interference (heart rate variability, car- 
diovascular autoregulation) is around 0.01 Hz. The respiratory 
fluctuation is around 0.2 Hz. The frequency band we chose for 
analysis (0.02-0.08 Hz) does not include these interferences and 
therefore, connectivity detection as a result of these interferences 
in unlikely. 

Bivariate methods in general can result in non-existing spu- 
rious connections when there is a propagation of information 
from one channel to others. A pairwise measure of connectivity 
will result in detection of connection between all possible con- 
nection pairs, ie. direct or indirect. However, since we are only 
looking to find channels which belong to the same network, use 
of a bivariate measure of connectivity can be justified. Most meth- 
ods for functional connectivity mapping in the literature based 
on fMRI or fNIRS also use seed based methods which is relying 
on the bivariate concept of finding coherence/correlation between 
channels with the seed channel (Greicius et al., 2004; Zhang and 
Raichle, 2010). 

In summary, the results of this paper suggest that the proposed 
method can be used to reveal underlying connectivity patterns of 
cognitive functions in the resting state through phase relations 
between hemodynamic changes in different brain regions. The 
results also indicate a left lateralization in the detected net- 
work which suggests the language system may be left lateralized 
already in newborns. Most fNIRS studies have shown left lateral- 
ization of the language network in older infants and in adults. 
However, data have revealed inconsistencies regarding lateral- 
ization of this network in newborn infants, when examining 
activation in response to language stimulation (Pena et al., 2003; 
May et al, 2011; Perani et al, 2011; Sato et al, 2011). The 
results of this study in terms of the lateralization of the network 



provide further information and evidence for characterization 
of the language network in newborn infants. Thus, these results 
significantly add to the theory on the neural foundations of 
language acquisition by revealing that language processing net- 
works are left lateralized even before language acquisition proper 
has begun. 
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